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^ ; Abstract 

«— ■«. . We consider a simplified model of a two-phase flow through a heterogeneous porous 

f***^ ' medium, in which the convection is neglected. This leads to a nonlinear degenerate 

parabolic problem in a domain shared in an arbitrary finite number of homogeneous porous 
media. We introduce a new way to connect capillary pressures on the interfaces between 
the homogeneous domains, which leads to a general notion of solution. We then compare 
/\i ' this notion of solution with an existing one, showing that it allows to deal with a larger 

class of problems. We prove the existence of such a solution in a general case, then we 
prove the existence and the uniqueness of a regular solution in the one-dimensional case 
for regular enough initial data. 

Keywords, flows in porous media, capillarity, nonlinear PDE of parabolic type. 

1 Presentation of the problem 

The models of immiscible two-phase flows are widely used in petroleum engineering, partic- 
ularly in basin modeling, whose aim can be the prediction of the migration of hydrocarbon 
components at geological time scale in a sedimentary basin. 

The heterogeneousness of the porous medium leads to the phenomena of oil-trapping 
and oil-expulsion, which is modeled with discontinuous capillary pressures between the 



different geological layers. 

The physical principles models and the mathematical models can be found in [?J [5J IT21 
1271 113) . The phenomenon of capillary trapping has been completed only in simplified cases 
(see [7]), and several numerical methods have been developed (see e.g. [T51I14I ). 

The aim of this paper is to introduce a new notion of weak solution, which allows us 
to deal with more general cases than those treated in P3], while it is equivalent to the 
notion of weak solution introduced in [14] on the already treated cases. We will consider 
a simplified model (|2l) defined page [SI in which the convection is neglected. 

We then give a uniqueness result in the one dimensional case which is inspired from 
the result in [7J and extends this latter one to more general situations, by requiring weaker 
assumptions on the solutions and applying to a larger class of initial data. 

We have to make some assumptions on the heterogeneous porous medium: 

Assumptions 1 (Geometrical assumptions) 

1. The heterogeneous porous medium is represented by a polygonal bounded connected 
domain il C K d with meas^d(n) > 0, where meosjRn is the Lebesgue's measure of 

W l . 

2. There exists a finite number N of polygonal connected subdomains (Qi)i<i<N of CI 
such that: 

(a) for all i G [1,N]. meas^d^i) > 0, 

N 

(b) \j a = n. 

(c) for {i,j) e [i, n] 2 with i jk j, a, n n 3 = 0. 

Each fli represents an homogeneous porous medium. One denotes, for all (i,j) G 
[1,N] 2 , r^j C 51 the interface between the geological layers fit and ilj, defined by 

Tij =dn i 'ndn j . 

We consider an incompressible and immiscible oil-water flow through SI, and thus through 
each H,i. Using Darcy's law, the conservation of oil and water phases is given for all 

(x,t)en i x{0,T), 

(j>id t Ui(x,t) - V ■ (r/ 0i i{ui(x,t))(Vp 0:I (x,t) - p g)) - 0, 

-4>idtUi{x,t) - V- (r)w,i(ui(x,t))(Vp Wt i(x,t) - pwg)) =0, (1) 

Po,i(x,t) -p w ,i(x,t) = Ki(Ui(x,t)), 

where m G [0,1] is the oil saturation in Qi (and therefore 1 — ui the water saturation), 
(pi G ]0, 1[ is the porosity of fli, which is supposed to be constant in each fti for the sake 
of simplicity, TTi(iii(x,t)) is the capillary pressure, and g is the gravity acceleration. The 
indices o and w respectively stand for the oil and the water phase. Thus, for a = o,w, p a ^ 
is the pressure of the phase er, r\ a _i is the mobility of the phase <r, and p a is the density of 
the phase a. 

We have now to make assumptions on the data to explicit the transmission conditions 
through the interfaces I\j: 

Assumptions 2 (Assumptions on the data) 

1. for allie [1,N], it, G C x ([0, 1],K), with n'^x) > for x e]0,l[, 




Figure 1: An example for the domain Q 



2. for all i G [1,N], r\ 0t i G C ([0, 1],R+) is an increasing function fulfilling n Oi i(0) = 0, 

3. for all i G [1,N], r) W) i G C ([0, 1],R_|_) is a decreasing function fulfilling n w ^{\) = 0, 
4-. the initial data uq belongs to L°°(Q), < uq < I. 

One denotes on = lim s _j.o ^i{s) and /3, = lim s _>.i 7Tj(s). We can now define the monotonous 
graphs iii by: 

( tt,(s) if s e]o,i[, 

7fj(s) = < ] - oo, ai\ if s = 0, (2) 

[ [p h +oo[ a a = i. 

As it is exposed in [14], the following conditions must be satisfied on the traces of Ui, p a ,i 
and ^Pa.i on T%j x (0,T), still denoted respectively Ui, p a ,i and Vp CT ,i (see [5]): 

1. for any a — o,w, (i,j) G [1,N] 2 such that T i: j ^ 0, the flux of the phase a through 



Tij must be continuous: 



r)<jA u i){^P<y,i ~ P°g) ' n * + r la;j{Uj)(^P<r,j - Pag) ' tlj = 0, 



(3) 



where n^ denotes the outward normal of IYj to $7* 



2. for any cr = o,w, (i,j) G [1, N] such that Tij ^ 0, either p CT is continuous or r] a = 0. 
Since the saturation is itself discontinuous across Tij, one must express the mobility 
at the upstream side of the interface. This gives 



Va,i( u i)(Pa,i ~Pa, 3 ) + ~ Va,j(Uj)(Pa,j - Pa,i) + = 0. 



(4) 



The conditions (j4j have direct consequences on the behaviour of the capillary pressures 
on both side of Tij. Indeed, if < Ui,Uj < 1, then the partial pressures p and p w 
have both to be continuous, and so we have the connection of the capillary pressures 
Tti(ui) = iTj(uj). If Ui = and < Uj < 1, then p / t > p 0t j and p w ^ = p w .j, thus 




71,(0) 



7L(1) 



Jti(l) 



1 

Figure 2: Graphs for the capillary pressures 



ftjiuj) < 7Tj(0). The same way, u, = 1 and < Uj < 1 implies ttj(uj) > 7^(1). If Ui = 0, 
Uj = 1, then p ,i > p j and p M .i < p w ,j, so 7T;(0) > 7Tj(l). Checking that the definition of 
the graphs -rti and ttj implies 7fi(0) (~1 7^(0) ^ 0, 7f » ( 1 ) n 7Tj(l) ^ 0, we can claim that (TJ| 
leads to: 

7Ti(«i) r\TTj(Uj) ^ 0. 



We introduce the global pressure in fli 



Vi{x,t) =p w ,i(x,t) 



Ui(x,t) 



Vo,i( a ) 



'o VoA a ) +Vw,i(a) 

(see e.g. [5] or [12]). and the global mobility in Qj 

»?o,i( , u»(a; ) *))»7iu,t(«t(a;,*)) 



7Tj'(a)da 



Ai(ui(x,i)) 



r]o,i{ui{x,t)) + ri w>i (ui(x,t)) 



(5) 



(6) 



(7) 



which verifies A-j(O) = A,(l) = 0, and Aj(s) > for < s < 1. Taking into account ([6]) and 
in ([T]), and adding the conservation laws leads to, for (x,t) € Oj x (0, T): 

<j>id t Ui{x,t) - V • (»7o,i(ui(a;,*))(YPi(:M) - p c g) - Ai(iti(»,t))V7Ti(ui(a;,t))) = 0, 

\er= o,w / 

(8) 

We neglect the convective effects, so that we focus on the mathematical modeling 

of flows with discontinuous capillary pressures, which seem to necessary to explain the 



phenomena of oil trapping. This simplification will allow us to neglect the coupling with the 
second equation of (|5]), and we get the simple degenerated parabolic equation in Oj x (0, T): 



fad t u t {x,t) - V- (Ai(u 4 (a;,t))V7r i («i(»,t))) =0 in ^ x (0,T). 



(9) 



In this simplified framework, the transmission condition ([3]) on the fluxes through Tij can 
be rewritten: 

X i (u i (x,t))V(iTi(ui(x,t)))-n i + \ j (u j (x,t)) , V(ir j (u j (x,t)))-n j =0 onTy x (0,T). (10) 

We suppose furthermore that Ui{x, 0) = uq(x) for x £ tti. In the remainder of this paper, 
we suppose to take a homogeneous Neumann boundary condition, The existence of a weak 
solution proven in section [3] can be extended to the case of non-homogeneous Dirichlet 
conditions. Nevertheless, homogeneous Neumann boundary conditions are needed to prove 
the theorem 14. 1[ and thus to prove the conclusion theorem 15.41 

Taking into account the equations ([3]), ©, (jTUJ) , the boundary condition, and the initial 
condition, we can write the problem we aim to solve this way: for all i S [1 , N] , for all 
j 6 [1,N] such that T^ ^ 0, 



(j>id t u t - V • {\i{ui)V-Ki(ui)) = 

7T,(Wi) r\TTj(Uj) ^ 

\i{ui)V{-Ki{ui)) ■ n, + Xj(uj)V(Trj(uj)) ■ n, 
Ai(uj)V(7ri(uj)) • n, = 
u,;(-,0) = U (x) 



inOi x (0,T), 
onT itj x(0,T), 
onr 4J x(0,T), (V) 

ondft.Ddn x (0,T), 
in O,-. 



Remark 1.1 All the results presented in this paper still hold if one not neglects the effect 
of the gravity and if one assumes that the global pressure is known, that is for problems of 
the type : 

(j>id t Ui + V • (qfi(ui) + \i{ui)(p - p w )g - Ai(Mi)V7ri(Mj)) =0 in Qj x (0,T), 
ni(ui) n nj(uj) y^9 on T it j x (0,T), 

^2 (q/fc( w fc) + X k {u k )(p - pw)s - A fc (ii fc )V7r fc (w fc )) ■ n k = on T tJ x (0,T), 

k=i,j 

(qfi(ui) + X l (u i ){p - p w )g - Xi(ui)V%i(ui)) ■ tii = on dVL t ndVlx (0,T), 

Ui(-,0) = u (x) in Q;, 

where fi is supposed to be a C 1 ([0, 1],K) -increasing function, Xi is also supposed to belong 
to C 1 ([0, 1],R+) and q satisfies 

• Vt, q€ (C 1 (n i x[0,T])) d , 

• V-q = inCliX (0,T), 

• q|ft. • rij + q| fii • n,- = on Tij x (0,T), 

• q • n = 0. 

In order to ensure the uniqueness result stated in theorem \5.1\ the technical condition (see 
Jjf or JMl) has to be fulfilled: 

V*, h o n\ Xi o ip- 1 g C°^/ 2 ([0, ^i(l)],R). 



Remark 1.2 In the modeling of two-phase flows, irreducible saturations are often taken 
into account. One can suppose that there exists Si and Si (0 < Si < Si < 1) such that 
Xi(s) = if s ^ (si,Si). In such a case, the problem (fP|) becomes strongly degenerated, 
but a convenient scaling eliminates this difficulty (at least if Si < Uo < Si a.e. in fli). 
Moreover, the dependance of the capillary pressure with regard to the saturation can be 
weak, at least for saturations not too close to or 1. Thus the effects of the capillarity are 
often neglected for the study of flows in homogeneous porous media, leading to the Buckley- 
Leverett equation (see e.g. ]19$). Looking for degeneracy o/iih) 7Ti(lt) is a more complex 
problem, particularly if the convection is not neglected as above. Suppose for example that 
TTi(u) = eu + Pi, where Pi are constants, and let e tend 0. Non-classical shocks can appear 
at the level of the interfaces I\- f j (see ]l(fj). Thus the notion of entropy solution used 
by Adimurthi, J. Jaffre, and G.D. Veerappa Gowda \Tj is not sufficient to deal with this 
problem. This difficulty has to be overcome to consider degenerate parabolic problem. But 
it seems clear that the notion of entropy solution developed by K.H. Karlsen, N.H. Risebro, 
J.D. Towers ]2(K \21[ \22j is not adapted to our problem. 

2 The notion of weak solution 

In this section, we introduce the notion of weak solution to the problem (fP|) . which is more 
general than the notion of weak solution given in |13[ [14] . Indeed, we are able to define 
such a solution even in the case of an arbitrary finite number of different homogeneous 
porous media. Furthermore, the notion of weak solution introduced in this paper is still 
available in cases where the one defined in P3] has no more sense. We finally show that 
the two notions of solution are equivalent in the case where the notion of weak solution 
in the sense of [2] is well defined. The existence of a weak solution to problem fT^J in a 
wider case is the aim of the section [3] 

One denotes by ipi the C 1 ([0, 1],R+) function which naturally appears in the prob- 
lem §P\ and which is defined by: Vs £ [0, 1], 

<Pi(s) = / A l (a)7r^(a)da. (11) 

Jo 

Remark 2.1 The assumptions on the data insure that tp[ > on ]0, 1[, and so we can 
define an increasing continuous function ip~ : [0, </?i(l)] — > [0, 1]. 

We are now able to define the notion of weak solution to the problem j73) • 

Definition 2.1 (weak solution to the problem (fP|) ) Under assumptions \7\ and[^ a 
function u is said to be a weak solution to the problem TP\) if it verifies: 

1. ue L°°(Q x (0,T))),0 < u < 1 a.e. infix (0,T), 

2. Vi £ [1,N], (fi(ui) £ L 2 (0, T; H (£li)), where m denotes the restriction of u to 

ax(o,r), 

3. TTi{ui)f\TTj{uj) / a.e. on T %] x (0,T), 



I for all^jeV{VL x [0,T)), 



"* iQi Jo 



N 



\ I I <piUi(x,t)dtip(x,t)dxdt + y / 4>iUo(x)i[>{x,0)dx 

j JHi JO 



(12) 



V/ / V^i(in(ar,t))-VV'(a;,i)(ia;dt = 0. 



The third point of the previous definition, which insures the connection in the graph 
sense of the capillary pressures on the interfaces between several porous media, is well 
defined. Indeed, since (pi(iii) belongs to L 2 (0,T; H l (ili)), it admits a trace still de- 
noted ifii(ui) on Yij x (0,T). Thanks to the remark |2~T1 we can define the trace of Uj 
onTy x(0,T). 

Remark 2.2 One can equivalently substitute the condition: 

ibis. TT%{ui) n TTj(uj) 7^ a.e. on Tij x (0,T), 
to the third point of the definition \2.1l where TTi is the monotonous graph given by: 

( mis) ifse]o,i[, 

*<(') = { ^ a ^°^ l/s = °' (13) 

[A,max(/3,)] if 8 = 1. 

We will now quickly show the equivalence between the notion of weak solution to the 
problem (fPJ) and the notion of weak solution given in |14| . in the case where this one is 
well defined, i.e. N = 2 and max(ai, a 2 ) = a < j3 = min(/3i, /3 2 ). We denote as in P3] the 
truncated capillary pressures by tti = max(a, 7Ti), tt 2 = min(/3,7r 2 ), and we introduce the 
problem %P\. which is treated in |14j . 

<p t d t u t - V • (Ai(ui)V7ri(ui)) = in O, x (0, T), 

7ri(wi) = 7r 2 (u 2 ) on T it j x (0,T), 

Ai(«i)V(7ri(«i))-nx+A2(«3)V(7r2(«a))-n2 = on T itj x (0, T), (7>) 

Xi(ui)V(Tn(ui)) • n, = on 30* n 30 x (0,T), 

«i(-,0) = u (x) in Oj. 

Then it is easy to check that: V(si, s 2 ) G [0, l] 2 , 

Ti-i (si) =k 2 {s 2 ) <^7ri(si)n7r 2 (s2) 7^ <=> tti(si) n 7f 2 (s 2 ) 7^ 0. (14) 

In order to recall the definition of weak solution, we have to introduce the function 

[a, f3] -^R 

* : i f P 

p n- / min (Xj o 7r (a)) da. 

J a J- 1 . 2 
Vf is increasing, and for i = 1, 2, ^ o 7^ o ^7" is a Lipschitz continuous function. 
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Figure 3: Truncated capillary pressures 



Definition 2.2 (weak solution to the problem (fP\l) A function u is said to be a weak 
solution to the problem \P\) if it verifies: 

1. u€ L°°(n x (0,T))),0 < u< 1 a.e. infix (0,T), 

2. V? G {l,2},^( Ui ) G ^(O^sif 1 ^)), 

5. w : fix (0,T) ->• M, defined for (x,t) G 0, x (0,T) byw(x,t) = *07r i (u i )(o;,t) belongs 
to ^{Q^T-H 1 ^)), 

4. for alliJjeV(Qx [0,T)), 



f T 
a Jo 



E 



JV 



cj>iUi(x,t)dttp(x,t)dxdt + 2^ / <t>iUo{x)ilj(x,0)dx 



N 

E 



Hi JO 



V(/?i(wi(cc, £)) • Vt/>(x, t)dxdt = 0. 



Remark 2.3 TTie notion of weak solution to the problem \P\) can be adapted in the case 
where there are N > 2 homogeneous domains, but we have to keep conditions of compati- 
bility on (ai)i<j< w and (/3,)i<j<jv- 

Proof of the equivalence of the weak solutions 

On the one hand, if u is a weak solution to the problem (lPJ) in the sense of definition 
12.21 then for a.e. t G (0,T), w(-,t) G i7 1 (il),and particularly w(-,t) admits a trace on r^-, 
whose value is in the same time ty(7ti(ui(-,t))) and $f(jtj(uj(-,t))). Since \& in increasing, 
for a.e (x,t) G T^j x (0,T), 7T,-fuj(a;, i)) = itj(uj(x,t)). Using (fl4|) . we conclude that any 
weak solution to the problem \P\ is a weak solution to the problem fP]) in the sense of 
definition 12.11 



On the other hand, if u is a weak solution to the problem §P\ in the sense of defini- 
tion 12.11 then thanks to (|14p . for almost every (x,t) £ Tij x (0,T), 

Ki(Ui(x,t)) = TTj(Uj(x,t)) ^\&07Tj O yj" 1 (^j(Uj(x, <))) = * O 7Tj- O (pj 1 fa ( U j (x, t))). (15) 

Since \I/ o 7Tj o y~ is a Lipschitz continuous function, the second point in definition 12.11 
insures us that $ o 7Tj(uj) belongs to L 2 (0,T, H 1 ^^) for i = 1,2, and (|15[) insures the 
connection of the traces on Ty x (0, T), then the third point of definition 12.21 is fulfilled 
and m is a weak solution to the problem \P} ■ B 

Remark 2.4 We can define a function tt~ ,i £ [1,N], which verifies tt~ o 7Tj(s) = s /or 
any s 6 [0, 1]. Using the function defined on M. by ^(p) = J_ min^i^Aj o jr~ (a))da, 
it is easy to check that we can equivalently substitute the function \P o 7Tj(mj) to Vt o 7Tj(itj) 
in the third point of dehnition \2. c A This function is still defined if a > /3, but it becomes 
identically 0. so the notion of weak solution to the problem \V\i is weaker than the notion 
of weak solution to the problem KP)) . Indeed, in such a case, u(x,t) = Uo(x) = a G]0,1[ 
for any (x, t) £ O X (0, T) is a weak solution to the problem \P\) . but it does not fulfill the 
third point in definition \2.1\ 

3 Existence of a weak solution 

The aim of this section is to prove the following theorem, which claims the existence of a 
weak solution to the problem (j^J. This result has already been proven in section [2] in the 
case N = 2 and a > /?, for which the notion of weak solution in the sense of definition ^. II 
is equivalent to the notion of weak solution in the sense of definition 



Theorem 3.1 (Existence of a weak solution) Under assumptions^ and\^ there ex- 
ists a weak solution to problem fP\l in the sense of definition^ 



Proof 

In order to prove the existence of a weak solution to the problem f23) in the sense of 
the definition 12.11 we build a sequence of solutions to approximated problems (|16[) , which 
converges, up to a subsequence, toward a weak solution to the problem (fPJ) . The approx- 
imated problems do not involve capillary barriers, so existence and uniqueness of such 
approximated solutions is given in [9] . We let the proof of the following technical lemma 
to the reader. 

Lemma 3.2 There exists sequences (Aj„)„, (^i, n )n belonging to (C°°([0, 1], IR)) N such 
that, for i £ [1, N], and for n large enough: 

• A i>«|[0,i/n]u[i-i/n,i] = ~2> X iA s ) > 2~~2> f or al1 s e I ) 1 ]; A v»i -> Aj uniformly on 
[0,1],' 

3 1 

• 7Tj, n (0) = 7Tj,„(0) -5- -00, 7Ti,„(l) = 7Tj in (l) -» +OO, K n 2 > n' i>n > -, 7Ti )n -> 7Tj in 

L 1 (0, 1), TTi iTl — \ iTi and n'^ n — > 71^ uniformly on any compact set o/]0, 1[, 

• i/ie function ipi >n : s 1— >L Ai. n (a)7r^ n (a)da furthermore fulfills tfi >n ([0, 1]) = </?j([0, 1]) 
anrf ^ )W — > ifi in W 1,o °(0, 1). 



We also define the increasing functions: 

( [a n> b n ] -+R 

The conditions on the functions on the intervals [0, —] U [1 — —, 1] insures that for any fixed 
large n, the functions (<p^ n °^n ° ^n 1 )' are Lipschitz continuous. Then thanks to [3J, for 
all n, the approximated problems: 

4>id t U t , n ~ V • {Xi,n{Ui,n)'V^i,n(Ui,n)) = in 0, X (0,T), 

^i,n(Ui,n) = ^j,n{Uj,n) On l\j X (0,T), 

A,>(«i,„)V(7rt 1 n(Mi,n))-lli+A ]7 -,n(W ]7 -,„)V(7rj>(Uj,„))'Ilj = On Tij X (0,T), 

A l ,„(Mi,™)V(7r ij „(u l ,„)) • n,- = on <9f^ n dfl x (0, T), 

u,, n (a;, 0) = u (a;) in fi. 

(16) 

admit a unique weak solution in the sense of definition 13.11 given below, and this solution 

belongs to C([0,T],ZP(fi)) for 1 < p < +oo. 

Definition 3.1 (Weak solutions for approximated problems) 

A function u n is said to be a weak solution to the problem hlb}) if it verifies: 

1. u n <E L°°(n x (0,T))),0 < u n < 1 a.e. infix (0,T), 

2. \/i e {1,2}, ^ n ( Ui , n ) G tf^T-H 1 ^)), 

3. w n : Q x (0, T) — > R, defined on £1; x (0, T) 6?/ «i„ = \]/ n o 7^ n (u, „) belongs to 
L\0,T-H\tt)), 

4. for all V' £»(!! x [0,T)), 

N „ „T N . 

y I I 4»iUi t n(x,t)dtip(x,t)dxdt +y^ / <f)iUo(x)ip(x,0)dx 
-V) / / V^i,n(iti,n(»>*)) • Vip{x,t)dxdt = 0. 



(17) 



= 1 



The proof of existence of a weak solution given in [5J, shows that for all i G [1, N], for all 
n, there exists Ci > not depending on n such that, for all i € [1, NJ: 

||Vi,n(Wi, n ) ||L2(o,T;jyi(fii)) < Cl ll^J.n I|li(0,1) ' ( 18 ) 

thus (yi, n (wi,n))n is a bounded sequence of L 2 (0, T; H 1 (fli)) using lemma I3~1Z1 A study of 
the proof of the time translate estimate used in [3J Q3], and detailed in [TCI lemma 4.6] 
leads to the existence of C2 not depending on n such that: 

\\(Pi,n(Ui, n {','+r)) ~ <Pi,n( u i,n('>'))\\L*(aiX(0,T-T)) ^ rC 2 |Ki,n||l,i(0,l) llv?LIU~(0,l)- ( 19 ) 

Using lemma f5TS1 once again, estimates (JT5J), (flU)) allow us to apply Kolmogorov's compact- 
ness criterion (see e.g. [8]), thus we can claim the relative compactness of the sequence 
(<Pi,n(ui, n ))n in L 2 ^ x (0,T)). There exists fc e L 2 (0,T;if x (r2i)) such that 

(fi,n(Ui,n) -> fi in i 2 (^ l x (0,T)), 

10 



<Pi,n(ui,n) ~> fi weakly in L^O,! 1 ;]? 1 ^)). 

Let us now recall a very useful lemma, classically called Minty trick, and introduced in 
this framework by Leray and Lions in the famous paper |24| . 

Lemma 3.3 (Minty trick) Let (<j) n ) n be a sequence of non- decreasing functions with for 
all n, <p n : R — > R, and let (f> : R — ► R fee a non- decreasing continuous function such that: 

• n — >■ </> pointwise, 

• i/iere exists g £ L ; 1 oc (R) swc/i i/ia£ |^ n | < g. 

Lei fee an open sufesei o/R fc , fc > 1. Lei (« n ) n € (L°°(0)) N , Zei u e L°°(0) and let 
f e L^O) suc/i i/iai: 

• u n ^f u in the L°°(0)-weak-k sense, 

L/ien 

/ - 0(«). 

Since < u.^ n < 1, (wi : „) n converges up to a subsequence to m in the L°°(Qj x (0, T))- 
weak-* sense. ((pi, n )n converges uniformly toward <p 2 ; on [0,1], and we can easily check, 
using Minty trick, that fi = (fi(ui) E L 2 (0,T; Lf 1 ^)). Thus we can pass to the limit in 
the formulation (|17[) to obtain the wanted weak formulation: 

AT T JV . 

2 / / 4 > iUi{x,t)dt'4){x,t)dxdt + 2_. I 4>i u o(x)ip(x,0)dx 
i=i J^i h i=1 Jcu 

N r pT 



^2 Vipi(ui(x,t)) ■ Vip(x,t)dxdt = 0. 

.._, Jo., Jo 



The last point needed to achieve the proof of theorem 13.11 is the convergence of the 
traces of the approximate solutions (wj, n ) ra on Tij x (0,T) toward the trace of Ui, and to 
verify that 7r,;(u,;) H ttj(uj) ^ a.e. on T^j x (0,T). 

Since fij has a Lipschitz boundary, there exists an operator P, continuous from iJ 1 (J7i) 
into H l (R d ), and also from L 2 (^) into L 2 (R d ), such that Pv\ n . = w for all u e L 2 (^). 
Then P is continuous from H s {VLi) into H s (R d ) for all s £ [0,1]. One has, for all v E 

H s (n t ), 



\Hh.m < \\Pv\\ H .oi<) < llP^m^PvV^^CMy^M^y 

One deduces from the previous inequality and from ([H?]) that for all s e]0,l[, for all 
r g]0, T[, there exists C3 not depending on n, r such that 

||<^,n(Wj,n(-, ■ +T)) - <£;,„(«;,„(-, •))ll|3(0,T-T;H'>(n i )) - ^"^ (20) 

For si > S2, Lf Sl is compactly imbedded in LP 2 , and then estimate (|2T)|) allows us to claim 
that the sequence (<Pi, n {ui,n))n is relatively compact in L 2 (0,T; H s (£li)) for all s s]0,l[. 
Particularly, one can extract a subsequence converging toward <Pi{ui) in L 2 (0,T; H s (Qi)). 
We can claim, using once again Minty trick, that the traces of ((Pi, n (ui,n))n on Tij also 
converge toward the trace of tpi(ui), still denoted <pi(ui) in L 2 (0,T;.H ra-1 ' 2 (I\ iJ -)), and 
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particularly for almost every (x, t) G I\.j x (0, T). Since (pi is increasing, (Mi, n (^ ) i))n 
converges almost everywhere on r^ x (0,T) toward Ui(x, t). 

Let us now check that 7fj(uj) fl ttj(v,j) ^ a.e. on Tjj x (0,T). For almost every 
(x,t) G Tjj- x (0,T) the sequence (TTi jn (ui,n(x,t))) n converges (up to a new extraction) 
toward ji(x,t) G M.. Since for all n, iri, n {ui,n(%it)) = ^j. n (uj in (x,t)), one has: 

ji(x,t) =Jj(x,t) a.e. on I\j x (0,T). (21) 

If Ui(x,t) G ]0,1[, then 7j(x, i) = 7r,(u,(a:, £)). If u,(a;,i) = 0, 7i(a;, t) < on, and 
7i(x,i) G 7fj(0). In the same way, if m(x,t) = 1, ji(x,t) G ffi(l). 

This achieves the proof of theorem 13.11 because relation pip insures the connection of 
the traces in the sense of: 

7Ti(wi) n 7Tj(«j) 7^ a.e. on r i;J - x (0,T). 



4 A regularity result 

In this section and in section [SJ we show the existence and the uniqueness of a solution 
with bounded flux to the problem fP|) in the one-dimensional case. We make the proofs 
in the case where there are only two sub-domains Hi =] — 1,0[ and VL2 =]0,1[, but a 
straightforward adaptation of them gives the same result for an arbitrary finite number 
of fij, each one with an arbitrary finite measure. We now state the main result of this 
section, which claims the existence of a solution with bounded spatial derivatives on Qi, 
where Q { = Q, t x (0,T). We also set Q =] - 1, l[x]0,T[ and T = {x = 0}. 

Theorem 4.1 (Existence of a bounded flux solution) Letuo G L°°( — 1, 1), < uq < 

1 such that: 

• ipi{u Q ) G W 1 ' 00 ^), 

• 7Tl (1*0,1) fl 7T2(U0,2) ^ on r. 

Then there exists a weak solution u to the problem TP\) such that d x tpi(ui) G L°°(Qi). 

All the section will be devoted to the proof of the theorem 14.11 As in section [3l we 
will get this existence result by taking the limit of a sequence of solutions to approximate 
problems (|16|) involving no capillary barriers, whose data fulfill the properties stated in 
lemma 13.21 

Proof 
We will now build a sequence of approximate initial data (uo,n) adapted to the sequence 
of approximate problems. 

Lemma 4.2 Let uq be chosen as in theorem ^. 1\ then there exists (uo. n )n such that, for 
all n, 

• < «o,n < 1, 

• 7Ti )n (lio,n,l) = 7T2,n{Uo, n ,2) On Y. 
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The sequence (uo,n)n furthermore fulfills: 

Km ||uo, n -Uo||oo =0, \\d x <Pi,n{v<),n)\\L'= o (ai)'^\\ d x<Pi('"'0)\\L< x '(n i )- ( 22 ) 

Proof 

Since 7Ti(uo,i) f~l 712(^0.2) 7^ 0, then there exists (ai in ,a2,n) £ [0,1] 2 such that one has 
7ri,n(oi,n) = 7i"2,„(a2,n) and |oi,„ - Uo,i| + l a 2,« - "0,2! -> 0. One sets, for x e Oj: 

U0,n(x) = ¥'„ (T Vi [ipi{Uo) + Vi,n{ai,n) - fi(uo,i)]) 

where 

f s if 86 [0,^(1)] = [0,^, n (l)], 

T -Pi( s ) — { fi.n0-) ^ s > fiO)> 

[ if s < 0. 

Then the sequence (ito,n) converges uniformly toward uq. For all n, < uo, n < 1 and 
either <9 x <£j,„(uo,n) = d x <pi(uo), or <9 x </?i,„(uo,n) = 0. ■ 

The approximate problem (|16[) admits a unique solution u n thanks to [5] , which belongs 
to C([0, T], L 1 (0)). Now, in order to get a L°°(Qi)-estimate on the sequence (d x ipi^ n (u n )) n , 
we introduce a new family of approximate problems (|23p for which the spatial dependence 
of the data is smooth. 

Let 9 e C°°(R),0 < 6> < 1, with 6{x) = if x < -1, and 6»(x) = 1 if x > 1. Let & e N*, 
one sets: 

• ^(z) = (1 - 6{kx))4>i + 0(kx)fa, 

• Xn, k (s,x) = - Q(kx))\l,n(s) + 6(kx)\ 2 ,n(s), 

• 7T„, fe (s,x) = (1 - 9(kx))iri >n (s) + 0(kx)7T2,n(s). 

We will now take a new approximation of the initial data. 

"o,n (fc=r (* + l)) if *<-!/*. 
«o,n,fc(a;) = S > , , ,.< 

Uo,n ( jt=T (»- £) ) ^ X>l/k. 

In the layer [— 1/fc, 1/fc], uo,n,fe is defined by the relation 

(1 - 0(fca;))7ri,, l (iio,n 1 A(aO) + ^(fcx)7r2,„(uo,r l ,A ; (x)) = 7ri,„(oi,„) = 7r 2 , ra (a 2 ,n), 

so that the approximate capillary pressure ^ n ,k(uo,n,k, •) is constant through the layer. 
Moreover one has either 

k 1 

\n,k( U 0,n,k,x)d x ( / jr n ,k( u O,n,k,x)) = _ <9 a ¥>i,n(^0,n) if |#| > 7i 

or 

9 x (7r„ ]fe (uo,n,fc,a;)) = if \x\ < -. 

So we directly deduce from the definition of Uo,n,fc the following lemma: 
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Lemma 4.3 Let n > 1, < Uo,« < 1 with ipi,n(uo, n ) 6 W^ 1,0o (fii) and 7Ti in (uo,n,i) = 
"■2^(10,11,2)1 i/ien i/iere exists a sequence (uo, n ,k)k satisfying, for all k > 2 ; £/W < 
"o,n,fc < 1 a«rf 

||A„,fc(U0, re ,fc,')9a;(Tn,fc(W0,n,fc, 0)l|oo < 2 max(||5 x ^,„(u ,n) ||oo), 

i— 1,2 

Wo,n,fc ~^ w 0,n * M ^ (fi) as k -^ +00. 

For any fixed fc > 2 and n large enough, we can now introduce the smooth non-degenerate 
parabolic problem (l2~3l) : 

cf> k (x)d t u n ^ k -d x {\„M(u n , k ,x)d x ir n . k (u nik ,x)) = 0, 

9xWn,fe(-l,*) = d x u n ,k(l,t) = 0, (23) 

u n ,fe(a:,0) = uo >nt k(x). 

Moreover, one can furthermore suppose, up to a new regularization, that wo,n.fc £ C°°([— 1, f ] 
Then ([2"B")| admits a unique strong solution u n , k G C°°([0,T] x [—1,1]) (see for instance 

HO [23]). 

Now one sets f n ,k(x, t) = \ n .k{u n .k, x)d x ir n>k (u ntkl x), so the main equation of (|23p can 
be rewritten: 

<f> k d t u n ,k = d x f n , k . 

A short calculation shows that f n .k(x,t) is the solution of the problem: 

9tfn,k = a n . k d xx f n . k + b n ,kd x fn,k, 

/n,k(-M) = /n,k(M) = 0, (24) 

fn,k(x,0) = A n] fe(u ,n,fe,-)^x(7I"n,fe(W0,n,fe,-)): 

where a n _ k ,b n>k are the regular functions defined below. 

a n ,k = K,k{Un,k,X) ttt— r > 0, 



,^[" n ,t(M n ,t ,a:) 

(/) fe (x) 



7 r\ Ml \ u x["n,k\ a n,ki -^ ) I , / nq 

Vfc = (A„,feJ (U n ,k,X) IfcZTi 1" ^n,k{ U n,k,X)C> x 



(^n,k)'(Un,k,x) 



4> k (x) 

The fact that wo.n.fc is supposed to be regular allows us to write the problem ([2~4")l in a 
strong sense (this is necessary, because this problem can not be written in a conservative 
form). In particular, f n ^ k satisfies the maximum principle, and thus 

||/n,fc||z,°°((-l,l)x(0,T)) < ||A„,fc(uo, n ,fc, O^^n.fc^O.n.fc; ')) IU=°(-1,1) • 

Thanks to the lemmas [4.31 and I4T21 we have a uniform bound on (f n ,k) : 

||/t.,*||l°°((-i 1 i)x(o,t)) < 2max(||a x ^(uo)||oo)- (25) 

Since the problem ([23)) is fully non degenerated (recall that A,;.„ > -^ and 7r,' n > — ) 
it follows that d x u n , k and dtu n . k are uniformly bounded respectively in L°°(Qi) and in 
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L 2 (0,T : H 1 (f2j)) with respect to k, then the sequence {u n ^)k converges toward u n in 
L 2 (Qi), and the limit u n fulfills, thank to estimate ([25]) : 

\\dxVi,n{u n )\\L<*>{Q i ) < 2max(||9a;Vi(uo)||oo)- (26) 

One has for all ^D([-l,l]x [0,T[), 

/ / cf> k u n , k d t ^ + / Ao^o - / / UM = 0. (27) 

Jo J-i J-i Jo J-i 

Thanks to $F5§, 

lim / / f n .kd x ip = 0. 
fc-H-oo,/ ,/_i 

One has u n .k — > u n in the i°°(Q)-weak- * and L 2 (Q) senses, uo.n,fc — > Mo,n in i 1 (— 1, 1) 
thanks to lemma [4731 Moreover, thanks to estimate (|25p . d x Tti^ n ^{u n ^) — > d x TTi in (u nt k) in 
the i°°(Q)-weak- * sense. Thus we can let k tend toward +oo in (|27|) to get 

ll2 < t > i v 'ndti>+ y] / 4>iU ,n1pO- V] / ^i,n(Un)d x Tti,n(Un)d x 1p = 0. (28) 

Vo l=lj2 Jn, i=l,2" / ^ ■ y ° i=l,2" /0 < 

Furthermore, using the fact that ir n ,k(u ny k,x) belongs to L 2 (0, T; i? 1 (Sl)) and, even more, 
that d x (ir n ^k(u n ,k, x)) is bounded uniformly in fc, we can claim that 7Ti. n (ui.„) = ^2,71(1*2, n)j 
and so u„ is the unique weak solution to the approximate problem (|16p for u 0n as initial 
data. 

When n tends toward +00, the sequence (u„) n converges, up to a subsequence toward 
a weak solution to the problem ^P\ , as seen in section [3J but the estimate (j^Hj) insures 
that 

d x <pi(u) e L°°(Gi). 



This achieves the proof of theorem 14.11 ■ 

5 A uniqueness result 

In this section, we give a uniqueness result in the one dimensional case in a framework 
where the existence results are stronger than the general existence result stated in the- 
orem [XT] Under a regularity assumption on the initial data uq, we proved in section 0] 
the existence of a solution having bounded flux, for which we give a uniqueness result in 
theorem 15. II and corollary [572] The bound on the flux will be necessary to prove that the 
contraction property is also available in the neighborhood of the interface {x = 0}. Then 
we show in theorem 15.41 the existence and uniqueness of the weak solution which is the 
limit of bounded flux solutions for any initial data uq with < uq < 1. Indeed, the set of 
initial data giving a bounded flux solution is dense in L°°(n) for the i 1 (f2) topology, and 
theorem 15. II has for consequence that the contraction property can be extended to a larger 
class of solution, defined for all initial data in L°°(fl). We unfortunately are not able to 
characterize them differently than by a limit of bounded flux solutions, and we can not 
either exhibit a weak solution which is not the limit of bounded flux solutions. 
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Theorem 5.1 (Z^-contraction principle for bounded flux solutions) Letu,v be two 
weak solutions to the problem FP]) for the initial data uq, vq. Then, if d x ifi(ui) and d x (fi(vi) 
belong to L°°(Qi), we have the following L 1 - contraction principle: Vi G [0,T], 

^ / (f>i(u(x,t)-v(x,t)) dx < ^2 / 4>% ( u o( x ) -v (x)) dx. (29) 

The first part of this section is devoted to the proof of the theorem 15.11 which, with 
theorem 14.11 admits the following straightforward consequence: 

Corollary 5.2 (Uniqueness of the bounded flux solution) For all u$ G L°°(— 1, 1) 
with < uo < 1, such that, for i = 1,2, <£>j(wo) G W 1,00 (^)> arl ^ ^i(wo,i) H ^2(^0,2) 7^ 0, 
i/sere exists a unique weak solution to the problem fPl) m i/ie sense of definition \2.1\ and 
such that d x ip t {u) G L°°(Qi); moreover u G C([0, T], L p (0)) /or all 1 < p < +00. 

Proof 

The proof of the theorem l5.1l is based on entropy inequalities, obtained through the method 
of doubling variables, first introduced by S. Kruzkov [18] for first order equations, and then 
adapted by J. Carrillo [11] for degenerate parabolic problems. Note that in the present 
setting, we only need doubling with respect to the time-variable, as it is done, for instance 
by F. Otto [55] for elliptic-parabolic problems (or in [5] for Stefan-type problems). 
In the sequel of the proof, we will only give the comparison 

2J / <j>i(u(x,t) -v(x,t)) + dx < 2J / <\>i ( u o(x) - v (x)) + dx. 

i=l,2 , ' Q i i=l,2 iJi 

The comparison with (•)"" instead of (-) + can be proven exactly the same way. 

Let u be a bounded flux solution to the one-dimensional problem, i.e d x ifi(u) G L°°(Qi), 
i = 1,2. The weak formulation of definition 12 . 1 1 adapted to the one-dimensional framework 
of the section can be rewritten, for all ip G V(Vt x [0,T[), 



' y, / <fnu{x,t)dt^p{x 1 t)dxdt + 2, I 4>iUo(x)ilj(x,0)dx 

-/ y2 d x (pi(u)(x,t)d x i/)(x,t)dxdt = 

Jo l=1 Js), 



(30) 



This formulation clearly implies, for % = 1, 2, for all ifj G C£°(fij x [0, T[) with ^(0, t) = 0, 

/ (f>iu(x,t)dtip(x,t)dxdt + / (piUo(x)'tp(x,0)dx 
Ja t t Jo* ( 31 ) 

/ d x <pi(u)(x,t)d x ip(x,t)dxdt = 
JCli 

Classical computations (see e.g. [SJ [TTJ [SB]) on equation (|3ip lead to the following en- 
tropy inequalities: for all weak solutions u,v, for initial data uq,vq, for all £ G T> + (£li x 
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[0, T[x [0, T[) such that £(0, t, s) = 0, 

r T i-T 



(f>i(u(x, t) — v(x, s)) + (dt£,(x, t, s) + d s £(x, t, s))dxdtds 
o Jo JCli 

t r 

/ (j>i(uo(x) — v(x, s)) + t;(x,0, s)dxds 

° T n ; (32) 

+< 



cf>i(u(x, t) — v$(x)) £,(x, t, 0)dxdt 

o JQi 

/ / d x (ifii(u)(x,t) - {pi(v)(x,s)) + d x £(x,t, s)dxdtds > 0. 
o Jo JUi 

Let us note here an important consequence of the entropy inequality (|32p (and of the 
corresponding one for (u — v)~ ), namely that u can be proved to satisfy 

ess — lim / \u(x. t) — un(x)\dx = . (33) 



Indeed, this follows by taking v as a constant in (|32[) and using an approximation argument, 
see e.g. Lemma 7.41 in [33] ■ We deduce the time continuity at t = for any solution and 
in particular for both u and v taken above. 

Now, let p £ C£°(R,R+) with supp(p) C [-1,1] and J R p(t)dt = 1. One denotes 
p m (t) = mp(mt). Let ^ e D+([— 1, 1] x [0,T[) with V(0, •) = 0. For m large enough, 
£(x,t, s) = ijj(x,t)p m (t — s) belongs to V + ([— 1,1] x [0, T[x [0,T[), and we can take it as 
test function in (j3"2")) . Then summing on i = 1 , 2 leads to 



T /-T 



o Jo , 

T 



y / <fii(u(x,t) — v(x, s)) + dtip(x,t)p m (t — s)dxdtds 



+ y^ / ^( u o(^) - v(x, s)) + ^(.T,0)p m (-s)dxds 

JO ,_, 2 JO, 

,r , (34) 



i=1.2 Jn 



X^ / 4> t {u(x,t) - v a (x)) + ip(x,t)p m (t)dxdt 



t ,t , 

I /] I 9 x ((fi(u)(x,t) - ifi(v)(x, s)) + d x ip(x,t)p m (t - s)dxdtds > 0. 

JO i=1 2 J^i 

We can now let m tend toward +oo in (|34|) . and using (|33[) for u and w, and the theorem 
of continuity in mean, we get: for all ip £ V + {Q, x [0, T[) such that i(>(0,t) = 0, 

/ >J / (f>i(u(x,t) — v(x,t)) + d t tp(x,t)dxdt 
Jo i=li2 Jai 

+ 5Z / &( u o( x ) ~ v o (x)) + ip(x,0)dx (35) 

»=1,2 "'^ 

^ / d x (tpi(u)(x,t) - ipi(v)(x,t)) + d x ip(x,t)dxdt > 0. 



We aim now to extend the inequality ([33)1 in the case where ^>(0, t) ^ 0, and particularly 
in the case i[i(x,i) = 6(t), so that the third term disappears in (|35l) . 
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To this purpose, let us set here Ui(t) = Mi(0, t) to denote the trace of Ui at the interface 
r (and correspondingly, Vi(t) = Vj(0,£)). We introduce the subsets of (0,T): 

• E u>v = {te [0,T] I ui(t) > vi(t) or u 2 (t) > v 2 (t)}, 

• E u < v = {t£ [0,T] | ui(t) < Vi(t) and u 2 (t) < v 2 (t)}, 
so that E u < v is the complement of E u>v in [0, T]. 

For all e > 0, one defines ip e {x) = max (l - ^,oY For all e P+([0,T[), we take 
(x,t) >-> 0(i)(l — ip E {x)) instead of ip(x,t) as test-function in ([55]) . thus we get: 



V / <j>i{u{x,t) -v(x,t)) + d t 8(t)(l - t/j e (x))dxdt 



dt>0. 



z—1 2 * 

0(4) f ( V i(«)(-e,t)-^(t;)(-e,t)) + -^i(«i)(t)-^i(«i)(*))' 
+( V2 («)(e,t) - ^(v)(e,t))+ - faj(u 2 )(i) - Mv2)(t))+ 



For almost every £ S E u < v , the function (ipi(u) — </>j(i>)) ('i i) admits a nil trace on {a; = 0}, 
thus the third term in the previous inequality can be reduced to the set E u>v obtaining 



i=1.2 



Y^ / <t>i{u{x,t) -v(x,t)) + d t 6(t)(l -ip e (x))dxdt 



n, 



+ Y. &(«o(s) - «o(aO) + (l - i>e)(x)6(0)dx 



Oi 



(36) 



i=1.2 



+ / 6{t) J^ / d x ((pi(u)(x,t) - ip l (v){x,tj) + d x ^ E (x)dxdt > 0. 



We show now the crucial point of the uniqueness proof, which is the subject of the 
following lemma. 



Lemma 5.3 For all 9 G Z? + ([0,T[) ; ifu,v are both bounded flux solutions, i.e. if one has 
d x ipi(u),d x (pi(v) e L aa {Q i ) one has, 

limsup / 0(i) V] / d x (<pi(u)(x,t) - ipi(v){x,t)) + d x ip s {x)dxdb < 0. 

e^O JE U>V !=12 JS!j 



Using the weak formulation (|30[) . we can claim that for any regular function fl S 

P([0,T[), 

lim / 0(t) V / d x (ipi(u) - ipi(v))d x ip e (x)dxdt = 0. (37) 

Since for i = 1,2, d x (ifi(u) — <fi(v)) belongs to L°°(fli x (0,T)), one has 



$(t) ^ / d x ((fi(u) - ipi(v))d x ip £ (x)dxdt 



<c\\n LH0 ,T), 



then a density argument allows us to claim that (J3"T|) still holds for any i? G L^O,! 1 ), and 
particularly for ■#(£) = #(£)llE u>v (t). Thus there exists A(e) tending to as e tends to 
such that 

6(t)J2 [ d x ((pi(u)(x,t)-(pi(v)(x,t))d x ip e (x)dxdt = A(e). (38) 

Splitting up the positive and negative parts of (ifi(u)(x, t) — (pi(v)(x, £)), (|38p becomes: 

9{t) ^ / d x (tpi(u)(x,t) - (pi(v)(x,t)) + d x ?p e (x)dxdt 

* =1 V 2 (39) 

0(*) 5~] / Q„(V4(tt)(K,t)-¥>i(t;)(x,t))-ax^e(aj)da;<ft + A(e). 

It is at this point that we actually use the monotony of the transmission condition, 
i.e. condition 3 in Definition 12.11 Indeed, the conditions TTi(ui(t)) fl ^(^(t)) ^ and 
7Ti(vi(t)) r\TT2(v2(t)) 7^ insure that : 

u\ > vi => U2 > V2 and u\ < v\ => U2 < t>2 • (40) 

Therefore, recalling the definition of the set E u>v and of ip e , the first term in the right 
member of (|39|) is non-positive, and then we conclude 

lira sup / 9{t) V] / d x (ipi(u)(x,t) - tpi(v)(x,t)) + d x ip £ (x)dxdt < 0. 

e->0 •/#„>„ i^l.2- 7 ^ 

This achieves the proof of lemma I5~3l and allows us to take the limit in inequality (f31>f for 
e -4 0. Then for all V G X>+([0,T[), one gets 

^ / (j)i{u{x,t)-v{x,t)) + d t ^{t)dxdt< Y^ / 4 (wo(x)-u o (a;)) + V'(O)dx. (41) 

One can also prove exactly the same way that 

J2 I (j>i{u{x,t)-v{x,t))-d t ^{t)dxdt< Y, I MMx)-vo(x))-ip(0)dx. (42) 

These inequalities still hold for ip = (T — t), and then if uq = v$, one has u = v almost 
everywhere in Q. Moreover we can take ijj{t) = U[ 0]S ](t) as test function in (|41[) to get the 
i^-contraction principle (|2T))) stated in theorem 15. II ■ 

In the sequel, we prove that for any uo in L°°(— 1; 1), < Uo < 1, there exists a unique 
weak solution of problem fP| which is the limit of a sequence of bounded flux solutions 
(«„)„, i.e. for all n > 1, d x ipi(u n ) G L°°(Q 4 ). 

Theorem 5.4 (Existence and uniqueness of the SOLA) Lei uo G L°°(— 1,1), < 
wo < I; <™^ iei (uo.n)n>i be a sequence of bounded flux initial data, i.e. for all n>\, 

• < U ,n < 1, 
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• <pi(uo, n ) e W^iQi), 

• 7Tl(«0,n,l) n7T 2 (wo,n,2) 7^ 

sitc/i i/iai 



lim \\uo n -u \\ L un) = 0. 

Let (u n ) n >i be the sequence of the bounded flux solutions to the problem fP\) for uo. n 
as initial data. Then the sequence (u n )n>i converges toward u in C(]0, T[, L p {— 1, 1)), 
1 < p < +oo ; where u is a solution to the problem KP\l , called Solution Obtained as Limit 
of Approximation (SOLA). Furthermore, ifu,v are two SOLAs, for initial data uq,vq, one 
has the following L 1 -contraction principle: Vi G [0,T], 

N . N 

V/ (/) l (u(x,t)~v(x 1 t)) ± dx<y2 (l> i (uo(x)-Vo(x)) ± dx. (43) 

This particularly leads to the uniqueness of the SOLA. 

Proof 

Let (uo, n ) be a regular sequence of initial data converging toward uq in L (— 1, 1) - one 
take e.g. uo tn G C£°(] — 1, 0[U]0, 1[). Then (uo, n ) is a Cauchy sequence, and thanks to (E!)l) . 
for alH G [0,T], 

AT . JV . 

y^ / 4>i\u„(x,t) -u m {x,t)\dx < y^ / &|uo,n(aO -wo,m(»)|da;. 

*=1 "fa* i=l "to* 

Thus («„)„ is a Cauchy sequence in C([0,T]; L 1 (f2)) and converges to a function u in 
C([0, T];!, 1 ^)). Since (u„)„ is bounded in L°°(Q), one has i^-m in C([0, T]; £ p (-l, 1)). 
We now have to check that u is a weak solution to the problem (|73) . It is easy to check, 
using to the L°°-bound of u n , that (fi(u n ) tends toward fi(u) in L p (ili x (0,T)), for all 
p G [l,+oo[. Thanks to (JTSJ), the sequence {<f>i{u n )) n is bounded in L 2 (0, T; iJ 1 (r2,)), and 
thus <pi(un) — > <Pi(u) weakly in i 2 (0,T;iJ 1 (il i )), and</3j(u„) converges in L 2 (0,T; H s (fli)), 
for all s G]0,1[, still toward <fii(u). Particularly, u n ^(t) tends toward Ui(t). Since the set 
{(a, b) G [0, l] 2 | 7Ti(a) fl ^(b) ^ 0} is closed, we can claim that 

7fi(ui0))n7r 2 (ii 2 (») ^0 for a.e. t€ [0,T]. 

We can also pass to the limit in the weak formulation in order to conclude that u is a weak 
solution to the problem (jPj, achieving this way the existence of a SOLA u. 

Let now v be another SOLA, obtained through a sequence (fo,n)n of regular initial 
data converging toward vq. Thanks to (|125)) . one has, 



y2 4>i\u n {x,t) -v n (x,t)\dx < y^ / l ?M'Wo,n(a:)-fo,n(a;)|d£, 

whose limit as n tends toward +oo gives the attempted .^-contraction principle: 

JV „ N 

/ / 4>i\u(x,t) — v(x,t)\dx < y / <fii\uo(x) — vo(x)\dx, 

4=1 ^ 1=1 •'«< 

and so the uniqueness of the SOLA, completing the proof of theorem 15.41 
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